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Abstract 

We study the evolution of fermion mass matrices considering the hypothesis of 
approximate flavor symmetries (AFS) in the Standard Model and a two-Higgs- 
doublet model. We find that the hierarchical structure is not significantly 
altered by the running, hence the assumption of AFS is entirely compatible 
with a grand unification scenario. 
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The presence of strong hierarchies among fermion masses and mixing angles is in general 
regarded as a signal of new physics beyond the Standard Model (SM). At present, many 
different models have been developed, searching for a mechanism which could give rise to 
the observed parameter hierarchies in a "natural" way. 

In this work, we will concentrate on a model based on the existence of approximate flavor 
symmetries (AFS), which has been recently proposed by Antaramian, Hall and Rasin [|TJ. 
This model assumes the presence of a set of global U(l) symmetries (one for each quark and 
lepton flavor) which are only slightly broken, conducing to different suppression factors in 
the fermion mass matrices. The possible extension to many Higgs doublets 0, as well as 
the presence of neutrino oscillations 0, have also been studied within this scheme. 

Our purpose is to test the possibility of inserting the idea of AFS into the context of a 
grand unified theory (GUT), where the SU (3)c ® SU (2) l<S>U (l)y gauge structure of the SM 
is nothing but the low-energy manifestation of a theory containing a single gauge coupling 
constant. In such a framework, it is natural to expect that the fermion mass matrix textures 
will be determined at the scale of grand unification (typically of order ~ 10 15 GeV), rather 
than the electroweak one. That is what would happen if the gauge and the U(l) flavor 
symmetries were spontaneously broken at the same time. The form of the mass matrices 
at the electroweak scale will then be obtained by considering the evolution of the Yukawa 
couplings through the renormalization group equations (RGE). 

Evolution of AFS with the energy scale. In the Standard Model, the Yukawa couplings 
can be written as 

C Y = £ (Ag QM-j* + Ag Q t D 3 ^= + Ag LiE^ + h.c. (1) 

where Qi and Lj are quark and lepton SU{2)l doublets and Uj, Dj and Ej stand for right 
fermion fields (we assume by now that only one Higgs doublet H is present). The hypothesis 
of AFS determines the order of magnitude for the coupling parameters A, leading to the 
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relations 



~ e 0t e^ , |Ag| « e Qi e . , |Ag| « e Li e £j (2) 

where the e^'s are dimensionless suppression factors originated with the flavor symmetry 
breakdown. Their values can be approximately obtained from the experimental data on 
quark masses and mixing angles. 

In refs. fj]-fl, the relations @ are assumed to hold at the scale Q 2 ~ M§. However, as 
stated above, the A matrices will evolve with energy. We are interested here in the form of 
the Yukawa couplings at the large energy scale Mx ~ 10 15 GeV, hence we will consider the 
one-loop RGE for the A's in the SM. These can be written as H 



A F (3) 



where 11 and Sp (as well as A F ) are matrices in flavor space, with F — U,D, E. We have 
used here the definitions 

t= 1 -HQ 2 /Ml) (4a) 

Tp = Tr (3 A a A ut + 3 X D X D ^ + \ E \ E] ) (4b) 

Su = -S D = \ U \ U ^ -\ D \ D \ S E = \ E \ E ^ (4c) 

Gu = 8g 2 + ^g 2 2 + ^g 2 , G D = 8 g\ + ~<? 2 2 + ^gf , G E = ^g 2 2 + ^g 2 (4d) 

being g l7 g 2 and g 3 the running coupling constants corresponding to the U(l)y, SU{2) L and 
SU(3)c symmetry groups respectively. As it is well known, their evolution is given by 

9*{h) = gKti) f 1 - #| (t2 - h) gfih)) 1 (5) 



57T 2 



with 



2D 1 4 43 4 

K 1 = — N+-, K 2 = -N-—, K a = -N- 11 (6) 
9 6 3 6 3 v; 



(AT = number of generations). 
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Let us write now the A matrices using their approximate form (|2|), and ignoring for 
simplicity the presence of complex phases. From (|3|), the e parameters corresponding to X u 
will evolve with the energy scale according to 

- 167r2 i ln ( e ^) =au ^ ( ? ) 

with 

au (t) = Gu(t) - T F (t) -§E4 fe £ (4, - 4.) (8) 

Similar expressions can be easily found for F = D and E. In this way, at the scale tx = 
l/21n(M|/M|) we find 

(eg,:^)^) = (e Oi et^)(0) exp (- fi* au(t) dtj (9) 

This means that the running of A^ (and equivalently X D,E ) has not changed the ratios be- 
tween the matrix elements: the sole effect has been the introduction of a global multiplicative 
constant. On the other hand, if AFS are assumed to hold at the large scale tx (i-e., just after 
the spontaneous symmetry breakdown), they will be preserved when the coupling constants 
run down to the electroweak energies. Thus the AFS hypothesis, as presented in ref. [0, 
should in principle be compatible with a grand unification scheme. 

It is worth to notice that the equation (0) is just an approximate result, since we have 
taken the relations (0) as if they were exact. In fact, this would conduce to rank 1 quark 
and lepton mass matrices, that is, to two zero mass eigenvalues for each fermion type. The 
corresponding (order unity) corrections, however, will appear only in the last term of au, 
being suppressed by the small e 2 parameters. 

Numerical analysis. Once the mass matrix elements are determined at the scale 
Q 2 = M|, the multiplicative constants (3f = exp(— l/(167r 2 ) Jq x apdt) can be estimated 
numerically. For definiteness, we will use the cqJcqj ratios given in ref. [[| (we assume that 
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they do not vary significantly when running from 1 to 90 GeV), together with the quark 
mass values at the scale of the Z mass [[|: 

e Ql /e Q2 ^0.2 e Q2 /e Q 3^0.04 (10a) 

m d ~5.6 MeV, m s ~ 0.11 GeV, m b ~ 3.4 GeV (10b) 

m u ~ 3.2 MeV , m c ~ 0.85 GeV (10c) 

We will concentrate in the quark sector (the lepton matrix X E approximately decouples in 
([/])), ignoring as before the presence of complex phases. The matrix elements \fj D can be 
approximately calculated from the values in (|T0"D taking into account the relations 

Tfl p. , . 

e .e fi « — ^ 11 

where ((ft) ~ 175 GeV is the vacuum expectation value (VEV) of the neutral Higgs field. 

Our numerical results are represented in figure 1, where the solid lines stand for the 
values of f3u and (3d as functions of the top quark mass. We also include a dashed line, 
which corresponds to the value of the matrix element Aj^ at the scale Mx- This parameter 
deserves special interest, since it is the largest of the A^ appearing in the Yukawa couplings. 
At the electroweak scale, the relation A^ 3 ~ m t /{(j}) conduces to A^ ~ 1 for the expected 
range of m t , hence the identification of eg 3 and ejj 3 as suppression factors is questionable. 
However, as it is shown in the figure, the situation is considerably improved at the Mx scale, 
where the value of Aj^ lies between 1/3 and 1/2 for m t varying from 150 to 200 GeV. To 
this respect, notice that although the symbol "~" in (|Tl"D indicates that the equalities are 
in general just approximate (within a factor 2 or 3), the large value of m t compared to m UjC 
implies that the relation has to be almost exact for F = U, i = 3. 



Two-Higgs- doublet case. Let us finally carry out a similar analysis for a two-Higgs-doublet 
model (THDM) with natural flavor conservation (NFC) |J. In this scheme, the form of the 
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Yukawa couplings will be entirely similar to that of eq. ([!]), except for the replacements 
H — ► Hi and H — > /f 2 - The NFC condition ensures that only one Higgs couples to all 
the quarks of a given electric charge, preventing in this way the presence of flavor-changing 
processes at the tree level[|. Notice that one has to deal with two non-zero Higgs VEVs, v i 
and v 2 , therefore it is usual to introduce a new parameter, 

tan/3 =— (12) 

Vl 

As before, we will assume the AFS structure (|2[) and study the evolution of the coupling 
"constants" with energy The relevant one-loop RGE for the THDM can be written as in 
(||) with the definitions [0] 

Su = X U X U ^ + h D X D ^ , S D = Va^ + X D \ DJf (13a) 
T v = Tr (3 X U X U ^) , T D = T E = Tr (3 A D A Dt + X E X E ^) (13b) 

being t, Se and Gf the same as in (|j). The evolution of the g^s has also the form (|5|), where 

the constants are now given by 

20 1 4 4 

K 1 = -N+-, K 2 = -N-7, Kz = -N-ll (14) 

In analogy with the Standard Model, it is immediate to find that the evolution of the 
A matrices does not change the relations between the matrix elements. We have for the 
THDM 

*S(**) = A S(°) ^ h) = A S(°) ex p (~i^3 jf^W dt ) ( 15 ) 

where 

= GuAt) - Tu, D (t) - \ £ 4 E ( 16 ) 



It has been argued that with the introduction of AFS, the NFC requirement could be avoided. 
However, in that case further assumptions would be necessary to suppress the appearing CP- 
violating phases ||. 
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Consequently, it is plausible to assume that the AFS are originated at a grand unification 
scale, just like in the SM case. 

In order to get numerical estimations for the pjjjj factors, we will take once again the 
values in fllUp. However, due to the two-doublet structure of the model, the relations (|TTJ) 
will now be replaced by 

e Qi e Di « — - e Qi e v . « — - (17) 

while the VEVs have to satisfy the constraint 

(yl+vff ~ 175 GeV (18) 

resulting from the electroweak symmetry breakdown. The value of tan (3 is a new unknown 
parameter, which could in principle be determined experimentally. 

In figure 2, we present our numerical results for ftjj (a similar behaviour is found for 
/3p h ^) and X^ 3 (Mx) as functions of m t , for different values of tan/3. It can be seen that if 
tan (3 ^$ 1, the value of is less than one only for m t ~ 150 GeV, growing severely as long 
as mt approaches 170 GeV. The regard of A^ 3 as a parameter reflecting an approximate flavor 
symmetry could be plausible only for large values of tan (3, where the Yukawa couplings for 
the up-type quarks are similar to those of the SM. Hence, taking into account the existing 
phenomenological bounds on tan/3 ||, the possibility of assuming AFS at the Mx scale 
within the THDM appears to be considerably reduced in comparison with the Standard 
Model case. 

I would like to acknowledge Prof. L. Epele for a critical revision of the manuscript. This 
work has been partially supported by Universidad Nacional de La Plata (Argentina). 
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FIG. 1. Numerical results for the (3jj,d factors (solid lines) and \^ 3 (Mx) (dashed) as functions 
of the top quark mass. 
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